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ABSTRACT 

Exact non-perturbative partition functions of coupling constants and external fields exhibit huge hidden symmetry, reflecting 
the possibility to change integration variables in the functional integral. In many cases this implies also some non-linear relations 
between correlation functions, typical for the tau-functions of integrable systems. To a variety of old examples, from matrix 
models to Seiberg-Witten theory and AdS/CFT correspondence, now adds the Chcrn-Simons theory of knot invariants. Some knot 
polynomials are already shown to combine into tau-functions, the search for entire set of relations is still in progress. It is already 
known, that generic knot polynomials fit into the set of Hurwitz partition functions - and this provides one more stimulus for 
studying this increasingly important class of deformations of the ordinary KP/Toda r-functions. 



1 Introduction 

For a long time integrability was thought to be an exotic phenomenon, useful for providing solvable examples 
in the textbooks and for enjoying the intellectual gourmands by mathematical beauties and possibilities to tie 
together the different branches of science. 

However, today, despite the growing importance of these motivations, especially the last one - actually, the 
key goal of the string theory [1], - integrability occupies a much bigger place in theoretical studies, and its role 
and significance is only increasing. The word "integrability" appeared hundreds of times only in the titles of 
the hep-th papers during last several years, well known subjects include matrix models, Seiberg-Witten theory 
of non-perturbative Yang-Mills fields, AdS/CFT-correspondence and many newly-emerging topics - Hurwitz 
and knot theories among the most recent ones. Ironically, when one performs this search in arXiv, it lists 
together the titles with "integrability" and "integrals", predominantly functional. There it is, probably, just 
the linguistics, but in fact this is symbolic, because integration is indeed the place, where the modern role of 
integrability has its origin. 

In fact, today we think [1] that integrability is the pertinent feature of non-perturbative (i.e. exact) partition 
functions - the quantities obtained by evaluation of functional integrals, and depending on the parameters of the 
integration measure (coupling constants), boundary conditions (background fields and vacuum condensates) and 
the integration domain (cut-offs). The main property of the integral is that it does not depend on the integration 
variables - and this is the property, which is finally reflected in integrability features of the partition function, 
often not immediately visible (hidden), when it is calculated by various techniques, which (most of them) do 
not explicitly respect integrability. 



2 Renormalization group 

Traditional way to think about integrable systems is in terms of commuting Hamiltonian flows on the phase 
space - but this is not at all straightforward to see and identify these flows in the study of non-perturbative 
physics. Perhaps, the most obvious place where at least one flow is studied for many years in this context, is 
the renormalization group. 

Namely, when one tries to define a partition function like 

Z(g k \A) = [ D<t>(x)e iS W H (1) 

J\p\<A 
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where, morally, 

S{4>} = J {{d4>f + Y,9k^)d d x (2) 

k 

it is possible to ask that the coupling constants g k depend on the cut-off A in such a way, that 

±- Z (g k {k)\k)=Q (3) 

This renormalization procedure is extremely interesting by itself, and itself involves a lot of group theory and 
integrable-like structures [2, 3] - but this is beyond the scope of this text. What matters, this principle defines 
a flow 

on the space of the coupling constants. Of course, it is just a single flow, to make it into a set of flows one 
should consider arbitrary deformations of the boundary (cut-off) of the integration domain [4, 5], what was 
never really studied in quantum field theory. 

Still, even looking at one flow, one can say something. The point is that the space of coupling constants 
is multi-dimensional, and in such case a generic flow is " chaotic" . This is what one could expect from generic 
renormalization group [6], but this is what was yet never observed in any QFT example. Partly this is reflected 
in various versions of Zamolodchikov's c-theorem [7], but actually the statement is stronger: the flow could be 
contractive, but still chaotic in orthogonal directions. The fact that this does not seem to happen, can be a 
signal of a hidden integrability - and in a context, much broader than any of the already well-studied examples. 



3 Ward identities and AMM/EO topological recursion 

Still, the flows are not the most convenient language to describe integrability in modern days. It usually shows 
up differently: as an infinite set of consistent differential equations on partition function Z(g k ). These equations 
are in fact nothing but Picard-Fuchs equations for the integrals, which depend on parameters - and they reflect 
the possibility to change integration variables, what can be alternatively imitated as the change of the coupling 
constants. In QFT such equations are also known as Ward identities. 

The archetypical example [8] is the matrix integral (matrix model) - with already many enough integration 
variables to observe non-trivial structures, and with still few enough to make as detailed calculations as one can 
wish. This calculation leads to the celebrated Virasoro constraints a la [9]: 

L n z(t) = Y,HT k + tk )—+h 2 Yl dt^r b =°> n>-i ( 5) 

k a+b=n 

which express invariance of 

Z(t k )= [ dM e i(-TrW(M)+-Z k t k Tr M *) (g) 
JNxN 

under the change 6M = eM n+1 . This allows to define the partition function as a D-module, for which (6) is just 
an integral representation. As to the integral, it can be well defined by introduction of background potential 
W{z) = — X^fc^fc 2 ^, usually polynomial, for example, quadratic. 

Amusingly, the dependence on this, originally auxiliary, potential appears to be one of the most interesting 
[10]. The terms of the genus expansion of the logarithm logZ of the partition function in powers of t'Hooft 
coupling constant h 2 N , are expressed as the meromorphic poly-differentials on a hyperelliptic Riemann surface, 

y 2 = W\zf f(z) (7) 

defined by the polynomial W(z) and some additional data f(z) - hidden in the naive definition (6). Virasoro 
constraints then turn into a recurrent set of relations, actually building a set of poly-differentials on the families of 
Riemann surfaces (often called spectral curves), which is now known under the name of AMM/EO topological 
recursion [10, 11] and which is now being found in a variety of different fields, from Seiberg-Witten theory 
of instanton sums [12, 13, 14] to knot theory [15]. From intcgrability-theory point of view this phenomenon - 
emergence of a new recursions from expansions of a Z)-module - is an example of emerging Whitham hierarchies, 
to be briefly described in s.6 below. 

In M-theory context the spectral curve can be nicely visualized as a surface, on which the 6d theory is 
compactified [23], and this picture leads to various deep insights [24], including the AGT conjecture [25]. 
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4 Group theory and quadratic relations 



Ward identities are linear relations, imposed on exact non-perturbative correlators. Topological recursion looks 
quadratic, but this is only a result of rewriting the action of quadratic operator on the exponential - the 
underlying equations are still essentially linear (though this is not so easy to see in some recent applications of 
the topological recursion). 

Remarkably, non-perturbative partition functions can also satisfy truly quadratic relations, which have noth- 
ing to do with linearity. They are usually called Hirota-like equations, and functions, that satisfy such equations, 
are known as r-functions. The best known is KP r-function, satisfying 



z k \ ( , z k 



dz e £*(**-t' fe )z-\ ^ fe + ±__ j T ^ _ t_ j = o (8) 

what, if expanded in powers of tk — t' k , provides an infinite set of quadratic differential equations (KP hierarchy): 
3(r 2 2T - t|) - 4(ti 3 t - tit 3 ) + (thht - 4rniTi + 3-r^) = 0, 

(9) 

The origin of such equations [16] is in the theory of Lie groups, where one can multiply representations 

Ri®R 2 = ® l R l (10) 

i.e. there is a comultiplication and intertwining operators. Quadratic relations appear when one consider two 
different decompositions with some common representation Q, belonging to both: 

R 1 0R 2 — >Q< — i? 3 ®^4 (11) 

Then a system of quadratic relations can be written on matrix elements of the group elements [17], in different 
representations - and this will be a Hirota-like equation for a generating function of such matrix elements. 
Unfortunately, not too many examples are worked out so far, beyond the KP/Toda family, which is associated 

to totally antisymmetric representations of the U(l) algebra. Remarkably, a new interest is now emerging to 
a quasiclassical limit of this construction for quantum groups - it is related to cluster algebras and to the 
Seiberg-Witten theory [18]. 

A very interesting fact about the tau-functions is that they can often be described in terms of auxiliary 
Riemann surfaces - spectral curves. In other words, the Liouville tori of integrable systems are often Jacobians 
of Riemann surfaces, and complex moduli are invariants of the Hamiltonian flows. It is still not clear, how 
general this phenomenon is, but it definitely holds for KP/Toda families [19]. 



5 Interplay between the linear and quadratic relations 

Thus non-perturbative partition functions naturally satisfy rich sets of linear Ward identities - as a result of 
the freedom to choose integration variables in the functional integral - and quadratic relations - as a result 
of existence of operator formalism in quantum mechanics, where amplitudes are matrix elements. With this 
understanding of their origins, it is clear, that the connection between linear and quadratic relations is not going 
to be very straightforward - as is the case of the subtle connection between the two formalisms in quantum 
theory. Indeed, even relation between the Virasoro constraints (5) and KP Hirota equation (8) for the simple 
matrix model (6) is not fully understood. 

In any case, it is clear that non-perturbative partition functions should be some kind of the tau-functions, 
but definitely not of generic type: in addition to quadratic they should satisfy at least one linear constraint 
- then all other Ward identities will follow. This single constraint is often called string equation, and such 
solutions to Hirota-like equations are called matrix-model t -functions [20]. 



6 Whitham hierarchies and Seiberg-Witten theory 

Normal in quantum field theory is not to define partition functions by integrating over all the fields completely, 
but rather do so with short-distance fluctuations, so that the answer remains dependent on the background 
slowly-varying fields, and effective action contains only terms with few derivatives. Such procedure gets espe- 
cially important, because we believe that the Standard Model of the fundamental interactions arises in exactly 
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this way from some still unknown UV-complete theory at Planckian scale (the superstring model [21] being one 
of the most popularized candidate, provided in the string theory context [1]). This procedure has, of course, a 
counterpart in integrability theory, where averaging over fast variables converts ordinary integrable hierarchies, 
satisfying Hirota equations, into effective theories of slow modes, described by the Whitham theory [22]. The 
crucial feature of Whitham hierarchy is that it depends on the choice of solution to original (high energy) hier- 
archy, and in the case when those solutions are parameterized by spectral curves,- it depends on the spectral 
curve. 

In other words, the belief is that the low-energy effective actions should be somehow described in terms 
of some additional hidden structure - auxiliary Riemann surfaces, or, better, the families of such surfaces, 
parameterized by the choice of the vacuum. Whitham dynamics defines the dependence of the low-energy 
effective action on the moduli - and the result is that it is often constrained by a new kind of non-linear 
equations: the WDVV-like system [26] 

TJj 1 T k = TkTj x Ti (12) 

where {Ti) . fc = d\- k logZ are matrices, made out of the third derivatives of the low-energy prepotential T = 
logZ, and derivatives are with respect to some special coordinates {a^} on the moduli space of spectral curves, 
introduced via the Seiberg-Wittcn (or special-geometry) equations 

I n = ai, I n = d l T= dl ° gZ (13) 

/a 4 M da, 

This construction appeared to be extremely successful [27] in application to the only known explicit example 
of exact non-perturbative effective action: the one obtained by summation over instantons in the N = 2 
supersymmetric Yang-Mills theory in 3,4,5 and 6 space time dimensions [13, 14]. In fact, it establishes a 
one-to-one correspondence between various quiver gauge theories and 1-dimcnsional integrable systems, like 
Lioville/Toda systems and spin chains, provides an effective description of S'-dualities [24, 28] and relates them 
to non-trivial spectral dualities between integrable systems, like [29, 30]. 

7 Lifting Whitham prepotentials to r-functions 

Of course, low-energy effective action is only a small remnant of original non-perturbative partition function. 
Likewise the Whitham prepotential is only a small remnant of original r-function, and there can be many 
different non-perturbative partition functions, leading to the same prepotential. It can be natural to look for 
the minimal among such extensions. 

In the case of the Seiberg-Witten theory, where Whitham dynamics is associated with entire world of 1- 
dimensional integrable systems, such an extension should be also distinguished. Today we already know that 
the corresponding minimal partition functions are closely related with conformal blocks of 2d conformal field 
theories - this is the celebrated AGT relation [25, 31]. The best possible formulation of it is in terms of Hubbard- 
Stratanovich duality in a peculiar Dotsenko-Fateev matrix model [32] , where Seiberg-Witten spectral curve and 
associated genus expansions are naturally appearing. A particular subset (particular Nekrasov-Shatashvili limit) 
of these conformal blocks is related to quantization of the Id integrable systems [33]. ^-duality relations seem 
to be realized in terms of 3d Chcrn-Simons theory. 

8 Chern- Simons theory and knot polynomials 

After distinguished Id and 2c? quantum field theories found their place in the world of universality classes of 
non-perturbative partition r-functions, the 3d Chern-Simons theory [34] is the natural next model to consider. 
And indeed the growing efforts are now applied in this direction. 

Chern-Simons theory is a version of Yang-Mills theory, with a topological action: 

f DAe U M M dA +^ 3 ) (14) 

and the observables are gauge invariant Wilson-loops 

TtrPcxpIa (15) 

J K. 
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along various lines (knots) K embedded into a three-dimensional manifold M.z. The Wilson- loop averages are 
often called knot functions, they depend on the space Ms, on the knot (or the link) /C, on the gauge group 
G, on its representation R and on the coupling constant h. For the simply connected Ms — R 3 or S 3 and for 
G = Sl(N) the knot functions are actually polynomials of the variables q = e 2lxlh and A = q N - and they are 
often called HOMFLY polynomials: 

H^f R f M *{q\A) = (j[Tr Rt Pexp£ Aj (16) 

There arc different ways to carefully define and evaluate HOMFLY polynomials, leading to various kinds 
of important quantities and structures, see [35] for classical results and [36] for relatively fresh reviews. In 
particular, if knots/links are represented by the closure of braids in the projection to the two-dimensional plane 
(this happens if the functional integral is evaluated in the temporal gauge Aq = 0), then knot polynomials 
exhibit a character decomposition [37, 38, 39, 40]: 

H%{q\A)=Y,C% Q {q) X Q{q\A) (17) 
Q 

which completely separates knot (K.) and group (A) dependencies. In the case of HOMFLY polynomials the 
expansion basis \Q is provided by Schur functions (the characters of Sl(oo)), and this implies a natural de- 
formation to MacDonald polynomials - such deformation of knot polynomials, depending on two parameters q 
and q' , indeed exists and leads [39, 42] to superpolynomials [41, 43]. Character expansion is actually defined on 
a bigger space of time variables instead of just A, where one introduces extended knot polynomials [39, 44], 

K(q, q'\tk) = £ C^ Q (q, q') XQ {q, q'\t k ) (18) 
Q 

which can be directly compared to r-functions. It turns out that - after appropriate summation over represen- 
tations R - they indeed give rise to r-functions in the case of the torus knots [44], but for generic knots this 
is not true: some modification is needed, cither of the extension of knot polynomials, or of the r-functions, 
perhaps, taking into account the existence of "quantum" parameter q. 

From the point of view of integrability, the most interesting seems the genus expansion near the point q = 1 . 
Exactly at this point the HOMFLY polynomial exhibits a very simple dependence on the representation: 

H%(q=\\A) = (o-£(A)) lRl (19) 

where \R\ is the number of boxes in the Young-diagram R, and □ denotes the fundamental representation, 
represented by a diagram with a single box. However, genus expansion is far more interesting [45]: the first 
correction is 

H%(q\A) = (o£(A)) W + {q-q-^ R {[2])^(A))^\f 2] {A) + ... (20) 

and in general 

H&q\A) = U(A)) exp (j> - q-^f 7%^U ( 21 ) 

\j \Q\<j+l M A )) / 

Here o~j\Q(A) are various special polynomials in A, with coefficients made out of Vassiliev invariants [46] of the 
knot K,. The i?-dependent coefficients (Pr(Q) are the characters of symmetric group Sao ~ which are nowadays 
studied in the framework of the Hurwitz theory. Note that for a given order j of the genus expansion only Young 
diagrams of the sizes |A| < j + 1 contribute, if one did not exponentiate the expansion, analogous restriction 
would be weaker: |A| < 2j. 

Extremely interesting seems the generalization of genus expansion (21) to superpolynomials, of which only 
the first-order term is conjectured [47, 48]: 

Pg(q, q'\A) = (P£{q, q'\A)) + [(q q' 1 ) v w (q' q'- 1 ) v R )of 2] {A) (^{A))^ 2 + ... (22) 

where </?_r([2]) = vr'—Vr, R' is transposed diagram R, and vr = J2jU ~^) r j- ^ should also include a non-trivial 
deformation of Hurwitz theory, describing refinement of symmetric group characters <pr(Q), which still remains 
to be found. 

An intriguing question is the consistency between (22) and the desired positivity of the colored superpolyno- 
mial: beyond non-(anti)symmctric representations the relation can be more involved [48], while (22) can instead 
describe MacDonald- inspired not- always-positive colored superpolynomials of [39, 42]. 
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9 Hurwitz partition functions 



The central object in Hurwitz theory [49] is the set of commuting operators W R , parameterized by Young 
diagrams R — {r\ > r 2 > ■ ■ ■ > 0} - labeling the conjugation classes of permutations of \R\ elements (see [50] 
for a non-commutative extension, associated with open rather than closed string theory). These operators have 
Sl(oo) characters (Schur functions) as their common eigenfunctions, while the corresponding eigenvalues are 
the symmetric group characters tpn(Q): 

WqXr = Vr(Q)xr (23) 
The simplest is representation of the operators in Miwa coordinates tk — ^TrX h , then 

^-■■u^( x ^y-- ( 24 ) 

In ordinary time- variables the first non-trivial is the cut- and- join operator 



W l2] =^(aMj b ^ b+{ a + b )t a+b ^) (25) 



which appears in various branches of science, including the construction of Khovanov knot polynomials [51]. 

Relation to Hurwitz numbers (counting ramified n-sheet genus-g coverings of Riemann surfaces) is through 
the Frobenius formula [52] 

Cover„, g (Qi, . . . ,Q P ) = ^ d 2 R 29 (p R (Q 1 )...Lp R (Q p ) (26) 

\R\=n 

It implies consideration of the generating functions of the form (at genus g = 0) 

Z(P\t) = exp PqWq^ e* 1 = ]T cxp ^ Q <p R (Q) j **(*) (27) 

which we call Hurwitz partition functions. 

It is clear from (21) that HOMFLY polynomials are actually the Hurwitz partition functions, evaluated at 
some knot-dependent point in the space of f3- variables. They can actually be rewritten in terms of the operators 
W, if one considers the Ooguri-Vafa generating function 

Z% v (t)=J2H R : XR(t) (28) 

R 

- note that this introduces a complementary set of times to (18). Then (21) implies that 

k 



Z§ v ® = exp fa, - q -V E W Q ) exp £ { - ~ > t k | (29) 

\j \Q\<j+i K(^)) / \k q q 

It was conjectured [53] that the character expansion of log Zov possesses non-trivial integrality properties (not 
to be mixed with integrability, at least immediately), which can be understood with the help of the theory of 
W operators [54]. 

Since operators Wq form a commutative algebra, Z((3\t) - and thus ZQ V (t), as its particular case - is clearly 
related to integrable systems. However, only for the unknot Z Q v turns out to be an ordinary KP r-function, 
satisfying Hirota equation (8) and possessing a free-fermion representation. For a general argument against 
such naive integrability for arbitrary knots (but not links) see footnote 3 in [44]. In fact, integrability properties 
of Hurwitz partition function are not fully understood and their connection to KP /Toda integrability remains 
obscure. In particular, KP integrability in t-variables is preserved only by the action of Casimir operators, i.e. 
when the Wq operators in (27) are linearly combined to form integrability-preserving Casimir operators, what 
does not happen for arbitrary values of (3q [55]. This can mean that either some change of t- variables is required 
to make integrability visible, or that KP-integrability is somehow deformed or generalized when one moves from 
2d conformal theories to 3d Chcrn-Simons. 



G 



10 Conclusion 



Intcgrability proved to be an extremely important part of quantum field theory: its main role is to describe 
full non-perturbative answers in quantum field theory - and this fact is deeply related to the very essence of 
quantum theory, whether it is formulated in operator or in functional integral formalism. It, however, remains an 
open question, what is exactly the right version of integrability theory which matters. Absolute majority of the 
known examples involves the simplest KP/Toda integrability, associated with the free fermion (antisymmetric) 
representations of the group U(l), though traces of more general quantum groups are already seen. Today 
the frontline seems to be at the study of Hurwitz partition functions, of which knot polynomials provide a 
vast set of examples: it is crucially important to understand and effectively describe their properties. Despite 
not immediately belonging to the set of KP/Toda r-functions, they seem to be at least their closest relatives, 
moreover, they seem to possess genus expansions, described in terms of the spectral curves and AMM/EO 
topological recursion. 
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